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Abstract 

A procedure for constructing bound state potentials is given. We show that, 

under the natural conditions imposed on a radial eigenvalue problem, the only 

special cases of the general central potential, which are exactly solvable and 

have infinite number of energy eigenvalues, are the Coulomb and harmonic 

oscillator potentials. 
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The celebrated Bertrand's theorem in classical mechanics states that the only central 
forces that result in closed orbits for all bound particles are the inverse square law and 
Hooke's law The extension of the Bertrand's theorem to spherical geometry led again 
to two potentials for which closed orbits exist 0. In the limit of vanishing curvature, these 
potentials reduced to the oscillator and Coulomb potentials. As is well known, these two 
potentials are also special in quantum mechanics. These are the only ones, which can be 
solved exactly in arbitrary dimensions and have infinite number of bound states. From 
the symmetry point of view, the degeneracy structure of these potentials admit a Lie group 
larger than the SO{d) group for a non-relativistic particle moving in rf-dimensional Euclidean 
space 0,11 . For the Coulomb problem, the group is S0{d + 1) for bound states and SO{d, 1) 
for scattering states. For harmonic oscillator, the corresponding group is SU{d). In case of 
spherical geometry, the eigenvalue problem can be solved exactly only for the two, above 
mentioned potentials [^,0 . These are the only central potentials for which the corresponding 
Schrodinger equations can be factorized to yield both the energy and angular momentum 
raising and lowering operators [[7|-[T(]|]. These two potentials are also special from the semi- 
classical point of view [|ri|]. Recently, it has been shown that, starting from a suitably gauged 
action, these two potentials can be derived in lower dimensions by appropriate gauge fixing 



I4 . However, unlike the classical case where a general central potential, under the twin 
constraints of bound and closed orbits, led to the Coulomb and oscillator potentials; in 
quantum mechanics, the natural conditions imposed on a radial eigenvalue problem, in a 
given dimension, have not yielded these two potentials as unique choices. 

In this paper, we construct a family of potentials, dependent on a parameter a, keeping in 
mind the two important features of the quantum mechanical bound state problems, i.e., the 
discreteness of their eigenspectra and the polynomial nature of the corresponding normaliz- 
able wavefunctions. Interestingly, under these general conditions, there are only two central 
potentials which are exactly solvable and have infinite number of bound states. These two 
unique potentials, appearing for the two values of a, one and two, are the Coulomb and the 
harmonic oscillator potentials, respectively. This quantum mechanical situation is an exact 



analogue of the Bertrand's theorem in classical mechanics. 

We begin with a Fock space spanned by the monomials, p"; where, n = 0,1,2,---. 
The most general diagonal operator which has a well defined action on this Fock space 
can be written as J^kL-oo^kD'^', here, D = and C^'s are constants. The two special 
cases corresponding to, (i) all the C^'s being zero except Ci = 1 and Cq = e (modulo 
an overall scaling factor) and (ii) all the C^'s being zero, except C2 = 1, Ci = (3, and 
Co = 6, can be unambiguously mapped to the Schrodinger equation by a series of similarity 
transformations. The first case leads to a class of potentials, dependent on a parameter a. 
It is explicitly shown that, the only two exactly solvable radial potentials, having infinite 
number of energy eigenvalues, are the harmonic and Coulomb potentials. Our method 
naturally gives their respective eigenfunctions. Further possibilities of solvable potentials 
are analyzed by performing a point canonical transformation. It is found that, no other 
exactly solvable central potential results and Morse potential arises as a conditionally exactly 
solvable case of this general potential. We then construct the second class of potentials and 
point out that, unlike the first, this case can not be solved as a radial eigenvalue problem, for 
infinite number of levels, for all the allowed values of a and for the various possible choices 
of the other parameters appearing in it. 

The choice (i), mentioned above, yields, 

(4 + .),W = ; (1) 

where, ri{p) = p~^ 

Performing a similarity transformation (ST): 0(p) = exp{yi/a;}?7(p), where, A = 
ap^~"^ + bp^'"-^ + cp~", a, a, b and c being arbitrary parameters, one gets 

^p^ + ap'-"^ + bp'-^^ + cp~" + e)(j) = . (2) 
dp dp^ dp 

The above choice of the ST is motivated from our desire to map (|1]) to the Schrodinger 
equation and also to have a closed form expression for the second order differential equation. 
Choosing = p""^^? it is straightforward to check that, 
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ap^ a ^ ^ dp 

+ -((e + a-2)p"-2 + (aa2 + (6-5a)a + 6a-26 + c)p-2)x = . (3) 
a 

For the purpose of clarity, we will cast as a radial eigenvalue equation in three dimensions; 
it can be easily seen that our method generalizes to (i-dimensions. Taking ^/^(p) = S{p)xip), 
one gets 

y dr^ r dr r"^ J ' 

where, the radial coordinate r = Ap, A being an appropriate length scale. Here, S{p) = 
p'^° exp{p"/ (2aa)}, Aq = (b — 6a + 2aa) / (2a) and / is the conventional angular momentum 
quantum number. 

The explicit form of the potential is 

V{r) =E + ^?ir2(--i) + g,r--' + g.r'' , (5) 

where, gi = ~gi\^^^~''\ g2 = h^^^"'', 93 = ~gz>?, ~9i = ny{8ma^X^), cj^ = Ti\2Aq - « + 5 - 
2e)/(4maA2) and ^3 = /i^{Ao(Ao + 1) - /(/ + 1) - [a(2 - a)(3 - a) + 6a - 26 + c]/a}/(2mA2). 

The corresponding unnormalized eigenfunctions, obtained through inverse similarity 
transformations, are 

V^(P) = p(^^-^)/^ exp {p7(2aa)} L^^/" (-p7(a«)) , (6) 

where, A = (1 — 6/a)^ — Ac/ a. Sign of should be chosen such that the resulting wave- 
functions are normalizable. The Laguerre polynomial, L^^/", is obtained by demanding 



that the (n+ l)th term in exp{y4/Q;}p ^" is equal to zero i.e., A"'~^^p = 0. This gives 



e^ = -an-(l/2)(l-6/a)±(l/2)VA . (7) 

Before proceeding to study the exact solvability of the special cases of the bound state 
problem obtained in (j^), it is worth repeating the corresponding scenario in classical me- 
chanics. In classical case, one demands the bounded motion of a test particle, subjected 



to a general central force of the form F{r) = —kt^ to trace a dosed orbit; this leads 
to two possible values for P, one and two [0]. These values give rise to the well known 
inverse- square and Hooke's laws, respectively. Interestingly, our general construction leads 
to a central potential, governed by a real parameter a. Much akin to the classical case, only 
for a = 1 and 2, as will be shown below, this potential can be solved for the complete set of 
energy eigenvalues. This is a quantum analogue of the classical Bertrand's theorem. 

The general central potential in (^) has free parameters gi or g2 for a = 1 or 2: one 
can then impose E + Qi = 0, for i = 1,2, in order that the potential is independent of 
constant terms and the corresponding Schrodinger equation is exactly solvable for all the 
energy eigenvalues. 

Case r. a = 1. 

To cast the eigenvalue problem in the standard radial form, i.e., to make = 0, we 
choose a = 1/(2(t), b = 2/a and c = [2 - l{l + l)]/{2a) with = -2m\^E/h^. We now 
demand that E + gi = (j^) then reduces to the Coulomb potential. 




Choosing g2 = — e^/(47reo), e„ in (0) gives the energy eigenvalues as 



Case 11: a = 2. 

Akin to the previous case, choosing a = l/(2cr), b = 1/a and c = — Z(Z + l)/(2cr) with 
a = —mujX'^/h, one can check that g^ = 0. Further demanding E + g2 = 0, the potential 
reduces to the harmonic oscillator potential: 



E = -gi = - 
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elh^ {n + l + 1)2 



V{r) = gir^ 



For, gi = muj^/2, the eigenspectra is 



E 



g2 = huj{2n + / + 3/2) 
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It is worth pointing out that only for these two values of a, one can make the potential 
independent of constants and hence obtain the full spectrum. We would like to add that, 



the above two cases can also be solved for gs |[T^-|T6|. For values of a different from 
one and two, the potential can be solved for E = 0, but for different values of the angular 
momentum quantum number 1. This situation is similar to the classical scenario, where, 
some values of (3 other than 1 and 2 may give rise to bounded motion, which is not closed. 

At this moment, one is naturally curious to see whether there exists the possibility of 
exactly solvable models emerging, under a point canonical transformation (PCT) p f{p) 
when a 7^ 1,2 |l^,[l8|. Performing the PCT, can be brought to the form (^, with a 
potential given by 
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(8) 



where, /' = df /dp. 

Different choices of / will give rise to different potentials; however, all of them may not 
be exactly solvable. For example the potential in (|) can be made independent of constants, 
even for continuous values of a, if one chooses /(p) = e^: 



2mX' 



= ^e'""' + —[b + a{a-l- 2e)]e' 
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(9) 



Since, the above potential contains the term '—/(/ + l)p^'^\ one can only solve the eigenvalue 
equation for a fixed 7'. This is the case of the conditionally exactly solvable potential 
encountered in Ref. |TB[. Physically this implies that, the potential is exactly solvable only 
in one dimension, z.e, without the centrifugal term. The corresponding energy eigenvalues 
are independent of /; they can be obtained from 



2mA2 
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and the unnormahzed eigenfunctions are 

One can easily see that, the Morse potential emerges as a special case for a = —1. 
Now, we proceed to find the second class of potentials starting from 



(l>' + /3L> + 5) 77(p) = . (10) 
Writing x{p) — 6xp{— 70}77(p), where O = cip"^ + bp°'~^, one gets 

{^'^''^'^^^'^^^Tp^^'^P^^^)^^^ ' ^^^^ 

where, F,(p) = Aip2(p"-i + A^f, F^ip) = Sip(p-i + B^^ + Bs and Fs(p) = Ci(p"-i + 
+ C3. The coefficients are given by Ai = A^^, = (a - 1)70, Bi = {aa + 2b)/{aAl), 
B2 = [2b + a(3 + /3 - 7)]/(2a^25i), B3 = 1 + (5 - B^Bl = b{aa - a + b)/{a^Al), 
C2 — b{(5 — 7 + 1) / {2aA2Ci) and C3 = — CiC|/6. Here, if one chooses A instead of O, 
then the resulting equation will contain still more higher derivative terms and can not be 
straightforwardly mapped to a Schrodinger equation. 

Now, writing x{p) — iFi{p)S{p))'~^ 'ijj{p); where, S{p) can be obtained from 

Ids _ F2ip)-2Fi{p) 1 
Sdp 2F,{p) p 

it can be checked that, ■0 obeys the Schrodinger equation with a potential, 

[5ip(p"-i + B2f + Bs][D,p{p--' + D2y + Ds] 



V{Xp) + 

4 



4Afp4(pa-l+^2)^ 

Di(p"-i + D2)[{2a - l)p-i + D2] 
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p2 + A2 J \p p"-^ + A2 ap"-! + ^2, 

[Ci(p"-i+>l2)' + C3 + (5] Z(Z + 1) 



Here, ^(p) = ^V{Xp), E = Di = (26 - 3aa)/{aAl), D2 = [2b + a(/? - 7 - 4a 

l)]/(2A2L'i) andDs^l + p- 4/L>i - D^Dl 



(12) 



Unlike the first class, these potentials can not be made independent of constants as a 
general radial problem, for any special values of a. Hence, the possibility of obtaining infinite 
number of levels, in arbitrary dimensions, does not arise here. However, it can be solved for 
the E = ^ state, for different values of /. 

At this point, we would like to add that, there have been many works in the literature, 
where point canonical transformation is used after writing 'i\}[x) = f{x)F[g{x)], for map- 
ping the Schrodinger equation to the confluent hypergeometric or hypergeometric equations 
||T9|^2|. Hence, from the known solutions of these equations, one is able to obtain solvable 
potentials in the Schrodinger eigenvalue problem, after imposing certain restrictions on f 
and g. In our case, we have started with the diagonal operators D = p-^ and Z)^ + aD 
acting on the space of monomials p". By inverse transformation, we, not only obtained the 
Schrodinger eigenvalue equation, but also the connection of the respective eigenf unctions 
with the monomials. In this approach, the special nature of the Coulomb and oscillator 
potentials, i.e., exact solvability with infinite number of levels in arbitrary dimensions, came 
out naturally without taking recourse to PCT. 

We would also like to point out that, a number of Calogero-Sutherland type interacting, 
many-body Hamiltonians, both in one and higher dimensions, can be mapped to the Euler 
operator J2xi'^ + c, through similarity transformations [^]. In fact, our partial motivation 



behind this work is due to a theorem established in Ref. 23 



All D dimensional particle Hamiltonians, which can be brought through a suitable 
transformation to the generalized form: H = J2iLi J2iLi xf^ can also be mapped to 

J2iLi J2iLi 3^F^^4o+^ exp{—d~^A}; where, the operator A is any homogeneous function of 
— 7JY and xf^ with degree d and e is a constant. For the normalizability of the wave functions, 
one needs to check that the action of exp{—d~^ A} on an appropriate linear combination of 
the eigenstates of YljLi YliLi x^P -^ij) yields a polynomial solution. It will be of great interest 
to apply the method employed here to the many-body systems and examine carefully, how 
far the results obtained for the single particle case generalizes to the N-particle systems. It is 
worth mentioning here that, apart from the Cologero-Sutherland model, which has recently 
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been mapped to harmonic oscillators the other many-particle exactly solvable system 
is of Coulombic type pl |. 
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